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Abstract 



Recent theoretical developments have shown that gauge mediation has a much 
larger parameter space of possible spectra and mixings than previously considered. 
Motivated by this, we explore the collider phenomenology of gauge mediation models 
where a general neutralino is the lightest MSSM superpartner (the NLSP), focusing 
on the potential reach from existing and future Tevatron searches. Promptly decaying 
general neutralino NLSPs can give rise to final states involving missing energy plus 
photons, Zs, Ws and/or Higgses. We survey the final states and determine those where 
the Tevatron should have the most sensitivity. We then estimate the reach of existing 
Tevatron searches in these final states and discuss new searches (or optimizations of 
existing ones) that should improve the reach. Finally we comment on the potential for 
discovery at the LHC. 



1 Introduction 



1.1 Motivation 

Supersymmetry is a promising solution to the gauge hierarchy problem. Its minimal realiza- 
tion - the MSSM - predicts a host of new particles with masses at the weak scale. These 
particles are presumably already being produced at the Tevatron, and they will be copiously 
produced at the LHC If they are someday observed in one of these experiments, then the 
questions of how SUSY is broken and how its breaking is mediated to the MSSM will become 
the central puzzles in particle physics. 

SUSY-breaking in the MSSM introduces ~ 100 new parameters beyond those of the 
Standard Model. The majority of these parameters introduce flavor-changing effects that 
are highly constrained by existing experiments. These imply that the mediation of SUSY- 
breaking to the MSSM must be essentially flavor blind. This is a puzzle in gravity mediation, 
where one must suppress flavor- changing effects from in principle arbitrary Planck-suppressed 
higher- dimension operators. On the other hand, it is an automatic consequence of gauge 
mediation [THO]. In this framework, the SUSY-breaking sector only communicates with the 
MSSM via the usual SM gauge interactions. Since these interactions are flavor blind, this 
guarantees that the parameters of the MSSM (in particular the squark and slepton masses) 
are diagonal in flavor space. 

Models of gauge mediation offer a rich variety of novel and distinctive collider signa- 
tures. Perhaps the most characteristic feature of gauge mediation is the light gravitino 
777.3/2 <C M wea k. This ensures that the effects of gauge mediation dominate over the Planck- 
suppressed contributions of gravity mediation. It implies that the lightest MSSM super- 
partner is actually the NLSP, which can always decay to the gravitino plus its SM partner. 
The gravitino will always escape the detector, leading to significant missing energy. Mean- 
while, the SM partner will tend to be central and energetic (assuming that the NLSP decays 
promptly), since it is typically much lighter than its NLSP parent. Given that pairs of NLSPs 
are produced in every SUSY event, it is clear from these considerations that the inclusive 
collider signatures of gauge mediation are largely determined by the nature of the NLSP. 

This fact has led to many searches for gauge mediation in the inclusive 77 + J£ T channel, 
especially after the discovery of anomalous Run I events at the Tevatron that contained 
ee77 + J£ T [7]. Searching for new physics in the inclusive 77 + J£ T channel was also well 
motivated from gauge mediation models of the mid 90's [HHTU]. In these "ordinary" or 
"minimal" gauge mediation models, there are weakly-coupled, ultraheavy "messenger" fields 
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$, $ transforming in complete multiplets of S77(5)gut which communicate with the SUSY- 
breaking sector at tree-level. In these models, the NLSP is typically either a bino or a 
stau, the former possibility giving rise to the 77 + J£ T signature. As we will review below, 
strong bounds have been set on these models by the CDF and D0 collaborations [TTI - 
[15] . and proposed searches in this channel have been put forth by the ATLAS and CMS 
collaborations [T5J[T7]. 

In recent years, there have been a number of theoretical developments in SUSY-breaking 
and model building which have greatly expanded the phenomenological possibilities of gauge 
mediation. For instance, in [T5] the most general renormalizable, weakly-coupled exten- 
sions of ordinary gauge mediation were analyzed in detail. It was shown that these "(ex- 
traordinary gauge mediation" models could realize very different spectra than those of 
ordinary gauge mediation. In one especially interesting set of models, the NLSP was a 
Higgsino-like neutralino, instead of bino-like as in ordinary gauge mediation. Apart from 
arising in this simple class of models, Higgsino NLSPs are interesting in their own right, as 
they are generic features of gauge mediation models which solve the little hierarchy problem 
|19j . That is, the neutralino NLSP is mostly Higgsino-like in the limit of small //, and this 
is also the limit in which minimizes the fine-tuning associated with the Higgs sector in the 
MSSM. 

In [2U] , the question of how to formulate gauge mediation in a model-independent way was 
examined. This resulted in the framework of General Gauge Mediation (GGM), in which 
the entire space of possible gauge mediation models can be investigated without having 
to rely upon specific models. The basic idea here was to view the SM gauge group as 
a weakly-gauged global symmetry of the SUSY-breaking hidden sector, and to relate the 
MSSM soft masses in gauge mediation to correlation functions of the associated conserved 
currents. The formulas obtained in this way are then exact in the hidden sector interactions 
but valid to leading order in the SM gauge couplings. Using this machinery, the complete 
set of independent parameters of gauge mediation can be derived, and in [20] they were 
shown to be: three arbitrary complex gaugino masses Mi^; and three real parameters 
controlling the five sfermion mass-squareds rug UD L E . It was also argued that in the GGM 
framework, the A-terms are always small. The extension to the Higgs sector was considered 
in [20] and in much more detail in [21]. Finally, in [22], simple toy models were constructed 
which demonstrate that the entire parameter space of GGM can be covered. This opens 
the door to new avenues of gauge mediation phenomenology. In particular, it implies that 
any MSSM superpartner can be the NLSP in gauge mediation. For some recent work on 
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the phenomenology of GGM, see [23H25] . Here let us briefly survey the possibilities. (For a 
similar survey that predates GGM and Tevatron Run II, see the last section of [26J.) 

If the NLSP is a colored sparticle, then the inclusive signature (assuming prompt decays) 
is dijets plus missing energy. Therefore, the bounds on this scenario are likely to be similar 
to those found in searches for MSUGRA, approximately 300 GeV [27] . From the point of 
view of naturalness, it is somewhat unappealing to have the bottom of the MSSM sparticle 
spectrum be so high. Also, this is an extremely well-studied final state at the Tevatron, so 
there is not much new that can be added. 

Next we come to slepton NLSPs. An important point to keep in mind is that there is very 
little direct production of sleptons at the Tevatron or LHC. Instead, sleptons are more likely 
to be produced from cascade decays originating in pair production of colored sparticles or 
charginos and neutralinos. So searches for sleptons will necessarily be more model dependent, 
involving more parameters of the spectrum. Furthermore, charged slepton NLSPs are already 
well-studied in the gauge mediation literature. Sneutrino NLSPs are an intriguing possibility 
which we will not discuss here; see however [28] . 

Finally, we arrive at neutralino and chargino NLSPs. Since these can be directly produced 
at the Tevatron, they are in some sense the minimal module for MSSM phenomenology at 
the Tevatron. However, charginos can only be NLSPs in narrow regions of parameter space, 
as we will review below, and as was analyzed recently in [29J. Thus we are left with general 
neutralino NLSPs as the remaining possibility for interesting, new, generic phenomenology. 

Motivated by the recent progress and these considerations, we will examine in this pa- 
per the phenomenology of promptly-decaying general neutralino NLSPs at the Tevatron. 
In general, the neutralino can be an arbitrary mixture of bino, wino and Higgsino gauge 
eigenstates. Therefore, having pairs of general neutralino NLSPs in every SUSY event will 
give rise to inclusive signatures involving any pairwise combination of 7 + J£ T , Z + J£ T or 
h + Fj t . While the 77+ J£ T channel has been extensively examined at the Tevatron, channels 
involving Z's and /i's have not been used to search for gauge mediation previously. (Prior 
to Tevatron Run II, the phenomenology of Higgsino-like NLSPs had been studied in detail 
|26 | l3"0Tf3"3"] . but not much experimental work has followed up on this.) Additionally, if there 
is a small splitting between the neutralino NLSP and the lightest chargino, there can be 
decay chains which end with W ± + J£ T which have been almost completely ignored in the 
gauge mediation literature. In this paper we will analyze these cases and others in detail to 
give a more complete picture. 

We choose to focus on the Tevatron and not the LHC in this paper, because the Tevatron 
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is a well-understood, smoothly-running machine which already has a large recorded dataset 
of about 6.5 fb _1 per experiment at the time of this writing (and more on the way). So it 
would be a shame not to push the existing Tevatron data as far as possible in the search for 
plausible new physics like gauge mediation. Also, using Tevatron data to constrain a large 
region of the parameter space of light neutralino NLSPs can help to focus future effort, both 
by model-builders and experimentalists. Of course, it is extremely important and interesting 
to ask how well the LHC can do, and how quickly it will surpass the capabilities of the 
Tevatron. These are issues that we will return to in a future publication [34] . 

In addition, the long-lived NLSP scenario is also clearly interesting. But since it involves 
qualitatively different physics and experimental techniques from the prompt case, we will 
examine it in a separate, forthcoming publication [35] . 

1.2 Search strategies and analysis methodology 

In the following sections, we will re-analyze existing Tevatron searches in final states relevant 
to general neutralino NLSPs and derive their constraints on the parameter space. We will 
also suggest, wherever possible, potential improvements to the existing searches that could 
enhance the sensitivity. Before moving on to the detailed discussion, let us first comment 
briefly on our general methodology. 

We will be focusing mostly on counting experiments. This is because SUSY cross sections 
at the Tevatron are generally rather low, so we are usually dealing with small numbers of 
events, making it difficult to study the detailed signal shape. Furthermore, we will be focused 
on obtaining current or projected 95% confidence-level exclusion limits, since these offer the 
most optimistic estimates for the "reach" of the collider. However, wherever possible, we will 
also mention our estimates of the 5a discovery potential and the 3a "evidence" threshold. 
We will use the maximum likelihood method to calculate the statistical significance of a 
given experimental search, following the discussion in [36]. For simplicity, we do not include 
systematic uncertainties in our estimates (in any event, we expect them to be small for 
signals involving leptons and photons). We also rely on tree-level estimates of signal rates. 
This is likely to lead to conservative estimates of the experimental sensitivity, because we are 
using accurate background rates, while the signal rates are likely to be enhanced significantly 
at NLO. 

When calculating signal rates, we will be assuming that both the colored sparticles and 
the sleptons are decoupled, so that the signal cross section comes entirely from pair produc- 
tion of charginos and neutralinos. Gluinos and squarks contribute significantly to the SUSY 
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cross section at the Tevatron only when m co i ore d < 450 GeV. Since it is often the case that 
they are much heavier than this in models of gauge mediation, it is reasonable to neglect 
them. Similarly, the cross sections for direct slepton production at the Tevatron are small 
unless the sleptons are extremely light, so even if sleptons were present in the spectrum, 
they would have little effect on our conclusions. These assumptions simplify the analysis 
considerably (fewer parameters to consider), and lead to more conservative bounds. 

All simulations discussed in the paper use Pythia |37J to generate signal and background 
samples, except for a few backgrounds with high-multiplicity final states (e.g. Zbb and bbbb 
in section [7]) for which we have used MadGraph [3B] to obtain the correct tree-level ma- 
trix elements (and, subsequently, Pythia to shower and hadronize the parton- level events). 
These samples were run through the PGS detector simulation [39], which includes models of 
tracking efficiency, calorimeter response, and 6-tagging. We have made small modifications 
of PGS to better account for the detector geometry and cracks at CDF and D0. Wherever 
possible, we have estimated backgrounds using PGS and compared them to published mea- 
surements and CDF or D0 background estimates that were computed with the full detector 
simulation and/or derived from data control samples. In channels involving real identified 
objects (leptons, photons) and real missing E?, we find that PGS is typically accurate at 
the 20% level, and is often much better. This gives us confidence that our signal simulations 
(which also involve real objects and real Fj t , of course) also achieve this level of accuracy. 
We find that PGS is less reliable when estimating backgrounds involving fake objects or fake 
missing Et', for these backgrounds, we will rely on published experimental results as much 
as possible. We further note that since the signal cross sections times branching ratios are 
generally rapidly falling functions of the sparticle masses, even a large error in the simulation 
of detector response can lead to a relatively small error in the estimated reach in mass. 

1.3 Relation to previous work 

In many respects, the work presented here revisits topics previously covered over a decade 
ago by many authors [261 130ti33] in preparation for Tevatron Run II. These important and 
comprehensive papers also studied the collider phenomenology of general neutralino NLSPs 
(along with many other cases), and they also used Monte Carlo simulations to examine the 
discovery potential in all the various final states. Since there is significant overlap with these 
older papers, it is worth highlighting briefly the ways in which our work here differs or goes 
beyond the previous literature. 

First and foremost, our reliance on detailed background estimates from both data and 
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Monte Carlo, and on the established capabilities of the Tevatron detectors in Run II, distin- 
guishes our analysis from earlier discussions of Higgsino NLSPs. Since these predated Run 
II, they were necessarily of a more rough and speculative nature. 

Second, we attempt to phrase things in a completely model independent fashion, in terms 
of the general neutralino parameter space (to be described later). This is in contrast with the 
old works which focused on a few "model lines" - one-dimensional slices through parameter 
space. 

Third, we also consider final states that were overlooked or discussed only briefly in the 
earlier literature; in particular, our analysis of wino co-NLSPs is an important addition. 

Finally, one of the main points of our paper is that - with the exception of 77 + Fj t - 
essentially none of the searches proposed in the old literature have been explicitly followed up 
on by experimentalists. At the same time, however, there have been searches in many of the 
important final states, either in a "signature-based" fashion, or motivated by scenarios other 
than gauge mediation. So in the present work, we hope to emphasize that experimentalists 
should be able to easily adapt existing searches to constrain an important portion of gauge 
mediation parameter space. Moreover, by optimizing these searches for general neutralino 
NLSPs, they may even have the possibility of discovery or 3a evidence. 

1.4 Outline 

The outline of the paper is as follows. In Section [2] we review the dominant decay modes 
of general neutralino NLSPs. We also examine in detail the special simplifying limits in 
which the neutralino becomes a pure gauge eigenstate - bino, wino or Higgsino NLSPs. In 
Section [3] we analyze the phenomenology of neutralino NLSPs at the Tevatron - production 
cross sections, the most promising final states (in terms of S/\fB), and the methodology we 
employ to analyze these channels. In section 4 we will discuss in detail the 77+ J£ T search and 
its applications to general neutralino NLSPs (including bino NLSPs, of course). In sections 
5-7, we will discuss searches relevant for wino, Z-rich Higgsino, and /i-rich Higgsino NLSPs, 
respectively. In Section 8 we summarize our results and conclude with some directions for 
future research, including comments about general neutralino NLSPs at the LHC. 
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2 Review of Neutralino NLSPs 



2.1 General neutralino NLSPs 

In this section we will review the different types of neutralino NLSPs and their decay modes. 
As discussed in the introduction, the general neutralino NLSP can be any linear combination 
of bino, wino, and Higgsino gauge eigenstates: 

4 

i=i 

where ip® = (B, W, H%, H®). The mass eigenvectors Nu depend on four MSSM parameters: 
Mi, the soft mass for the bino; M 2 , the soft mass for the wino; fi, the supersymmetric Higgs 
mass term; and tan (3, the ratio of the up-type to down-type Higgs VEVs. We assume the 
mass parameters are all real to avoid the SUSY CP problem; however, they can have either 
sign. (For more detailed discussion of the neutralino spectrum and parameter space, see e.g. 
Martin's excellent review of the MSSM [30].) 

The general neutralino NLSP has three possible decay modes: gravitino + {^/Z/H). As 
the parameters (Mi, M 2 , /i, tan/3) are changed, different decay modes will dominate. The 
general formulas for the decay widths are [3T| 132] : 



r(x?^G + 7 ) = \N u c w + N 12 s w \ 2 A 



1 , „ „ ,o\ I m 2 



r(x?^ G + Z) = ( \N 12 c w -N llSw \ 2 + -\N 13 cp-N liS p\ 2 ) ( 1 - ) A 



i 



2 / V x?. 

Tix^G + h) = \\N lz c p + N u sp\ 2 (l-^f\ A (2.2) 
1 \ m x ( 1 ) / 

These formulas are valid when the Higgs and Z decays are on-shell. (When necessary, we 
will also employ the full expressions including finite-width effects, which can be found in 
[33].) In the third equation, we have also assumed the SM Higgs decoupling limit. A is a 
dimensionful parameter that sets the overall scale of the NLSP lifetime: 

m l? / m*o \5/100TeVy' 



A 16vrF 2 ~ (lOO Gevi \ ^Fq ) 0.1mm' ^ 2 ' 3 ^ 
Here a/Fq is the fundamental scale of SUSY breaking; it is related to the gravitino mass via 
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m 3/2 = Fo/(y/SMpi). The range of possible y/To is roughly 

10 TeV < ^Fq < 10 6 TeV (2.4) 

where the lower bound comes from the requirement of a viable superpartner spectrum, and 
the upper bound comes from requiring gauge mediation to dominate over gravity mediation. 
This spans the range from prompt decays to lifetimes significantly longer than the size of 



the detector. Indeed, for much of the window (2.4), i.e. if y^o is somewhat large compared 
to 100 TeV, neutralino decays are macroscopically displaced. As discussed in section [TJ in 
this paper we will only focus on the case where the decay of the NLSP is prompt. 

Since the complete neutralino parameter space is four dimensional, it can be difficult 
to discuss the collider phenomenology in full generality. Instead, in this paper, we will 
mostly be studying various simplifying limits where the neutralino becomes a pure gauge 
eigenstate, i.e. where it becomes predominantly bino, wino, or Higgsino-like. For much of the 
parameter space, this is a pretty good approximation. Moreover, each limit has distinctive 
phenomenology. By considering each in turn, we capture the essential physics of general 
neutralino NLSPs. 



2.2 Bino NLSP 

The neutralino NLSP becomes bino-like in the limit \N\\\ ^> |-/V 12 |, | iVi3 1 , | -ZV14 1 . In terms 
of the MSSM parameters, this corresponds to \M\\ <C and \M\\ < |M 2 |. Substituting 



Nu -> 1 into ([2^2], we find 

r(x?^G + 7 ) = c 2 w A 

T(xi^G + Z) = siv^l-^fj A (2.5) 
r(x°-G + /i) = (2.6) 

with m^o pa \M\\. So we see that the bino-like NLSP decays to photons at least Cyy ~ 76% 








of the time, and to Z bosons the rest of the time. The branching fraction to photons can be 
further enhanced for lighter NLSPs, due to the m z phase space suppression factor. 
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2.3 Wino NLSP 



The neutralino NLSP becomes wino-like in the limit |ATi 2 | ^> |iVu|, \N\s\, \Nu\. This corre- 
sponds to |M 2 | <C fi and |M 2 | < \Mi\. In this regime, the decay widths are 

r(x?->G + 7 ) = s 2 w A 

Fix^G + Z) = ^(l-^f) A (2.7) 

V m x°J 

T(xi^G + h) = (2.8) 

with m^o « |M 2 |. Here the NLSP decays to photons at least 23% of the time, and to Zs the 
rest of the time. As with bino NLSPs, the branching fraction to photons can be enhanced 
if the NLSP is light, due to the Z-boson phase space factor. Indeed, for NLSP masses less 
than 150 GeV, the branching fraction to photons is greater than 50%. 

One special feature in the wino NLSP limit is that the lightest chargino - whose mass 
is also approximately |M 2 | here - is generically a co-NLSP. The decay rate of wino-like 
charginos to W + plus gravitino is given by 




T(xl^W + + G) = 1(1-^) A (2.9) 



xi 

If the chargino-neutralino mass splitting is sufficiently small, 

Am = |m-+ - m^>| < 1 GeV (2.10) 

then the decay length of xt ~^ Xi + ^* is longer than ~ 0.1 mm and the decay to gravitino 
(assumed here to be prompt) dominates [HHS]. When this is the case, we will refer to the 
wino-like chargino as a co-NLSP. 1 - 1 ) 

As is well known [IS], the splitting for wino charged and neutral mass eigenstates is 
abnormally small, due to a cancellation at leading order in mz/fi- The tree-level splitting 
goes like Am ~ m|//i 3 , so as long as /i is moderately large (/i > a few hundred GeV), the 
splittings will be sub-GeV. For this reason, when discussing the collider signatures of winos 
below, we will always assume that they are co-NLSPs. 

^We note that regardless of whether it is a co-NLSP, the lightest chargino will always decay promptly as 
long as the neutralino does as well. This should be contrasted with other scenarios with degenerate winos 
in the MSSM, for instance as in AMSB El. 
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2.4 Higgsino NLSP 

Higgsino NLSPs have been studied previously in [3014331, US]- Here we will review their 
basic features, with a greater emphasis towards the model-independent parameter space. 
The neutralino NLSP becomes Higgsino-like in the limit |iV 13 |, |iV 14 | 3> |iVu|, |iV 12 |, which 
corresponds to <C \Mi\, |M 2 |. In fact, this limit includes two distinct possibilities: 

N 13 = - V N U = -L (2.11) 

with 7] = ±lS 2 > The decay widths are given by 

r /~o n,\ l t , , 2 ( c w s w {M l - M 2 )m z \ 2 

r(x?-G + Z) = l(sp + vcpf[l-^f] A (2.12) 

r(*2-G + /i) = 

In the first equation, we have kept the leading-order term in the large Mi j2 limit, since this 
can affect the branching fractions at small NLSP mass where the Z and h decays are strongly 
phase-space suppressed. 

In figure [I] we plot the branching fractions for different values of tan j3 and 77. (For the sake 
of argument, we have taken m^ = 115 GeV in all these plots, and also for the remainder of 
the paper. Similar plots and discussion of the branching fractions can also be found in [30].) 
We see that for low mass NLSPs, there can be a significant branching fraction to photons, 
due to the phase space suppression factors (or simply being kinematically disallowed) in the 
other channels. Here we should emphasize that Br(7) at low mass depends sensitively on 



M 1 M 2 




4 








4 


m\ \ 


•4 



Mi and M 2 , through (2.12). For intermediate mass NLSPs, the decay is dominantly to Zs. 
Finally, for high mass NLSPs, the decay is still dominantly to Zs, except for rj = — 1 and 
small tan/3, where it is dominantly to Higgses. (For extremely high masses, one can also 
have Br(Z) ~ Br(/i), but unfortunately this regime is inaccessible at the Tevatron.) Given 
the different possible behaviors of the branching fractions, we will find it useful in this paper 



( 2 'ln terms of the underlying MSSM parameters, one finds 77 = sign(/i) x sign +tan 2 6*yi/j- While 

this depends on the relative sign between M\ and M2, in this paper we will only analyze the case that 
M\,M2 > 0. This is sufficient to capture the phenomenology of all neutralino NLSPs, and in this case, 77 
reduces to the more familiar variable sign(/i) used in many phenomenological analyses. 
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Figure 1: Branching ratios of the Higgsino NLSP to photons, Z's and Higg function of tunlsp, f° r 

T) = ±1 and tan/3 = 1.5, 20. In all the plots, M 1 = 500 GeV, M 2 = 1000 GeV, and m h o = 115 GeV. 



to distinguish between photon-rich, Z-rich and Higgs-rich Higgsino NLSPs. 

Similar to the wino NLSP, the lightest neutralinos and chargino Xii X2 an d Xi a ^ have 
masses ~ |/i| in the Higgsino NLSP regime. However, the difference here is that in practice 
one typically needs unrealistically large values of Mi, M 2 in order for them to be co-NLSPs. 
This is because (in contrast to the wino case) the splittings come at leading order in the 
small mz expansion, i.e. Am ~ m 2 z /Mi^- So in order for Am < 1 GeV, one needs Mi } 2 ^ 10 
TeV. For more realistic values of Mi and M 2 , the splittings are large enough that the heavier 
Higgsino-like states can decay promptly to the NLSP via off-shell W's and Z's going to jets 
or leptons. So we will not consider the possibility of Higgsino co-NLSPs in this paper. 

Another possibility that we will briefly mention is that of the chargino NLSP. This oc- 
curs in a special corner of the Higgsino limit, where /i < and — (p^J tyy < ¥r < 
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i i s , / ^w- The phenomenology of chargino NLSPs has recently been discussed in [2"5] . 
Since we are focusing on neutralino NLSPs in this paper, we will not discuss this scenario 
further. 



3 Overview of Tevatron Phenomenology 

In this section, we will give a broad overview of the phenomenology of neutralino NLSPs 
at the Tevatron, starting with their production cross sections. NLSPs are produced in 
every SUSY event, either directly or at the bottom of any decay chain. Thus the rate for 
NLSP production is just given by the total SUSY cross section times the relevant branching 
ratios. As discussed in the introduction, we are assuming that the SUSY cross section at the 
Tevatron is given entirely by the chargino/neutralino sector. Shown in fig. [2] are the total 
cross sections for bino, wino and Higgsino NLSPs, as well as the rates in various final states 
determined by the decay modes of the NLSP. Note that in the case of binos, the dominant 
cross section is not for direct production of binos themselves (this proceeds only through 
mixings), but for heavier Higgsino or wino states which then decay down to bino. The 
plotted cross sections and other estimates discussed in this section do not fold in geometric 
acceptances, except where otherwise noted. (Such acceptances will, of course, be included 
in the more detailed simulations in the remainder of the paper.) 

We clearly see from fig. [2] that regardless of the final state, the Tevatron will not be able 
to probe chargino masses above 300 GeV. (Of course, this depends on the assumption of 
decoupled colored sparticles.) Moreover, to probe this entire range one needs to have high 
signal acceptance and/or no branching fraction suppressions, as well as to be able to suppress 
the SM background entirely. As we will see below, only the bino-like NLSP has a chance of 
satisfying all these conditions, through the 77 + J£ T final state. 

Now let us look at the final states for the different types of neutralino NLSPs and discuss 
which are most relevant. In the case of bino NLSPs, it is clear that an inclusive 77 + p T 
signature is the best search channel. The branching fraction to photons is large, photons 
are clean experimental objects, and simple cuts in this channel can eliminate essentially all 
SM backgrounds. This channel is useful even when the bino content is highly mixed with 
Higgsino, as we will show in section |4j 

For wino and Higgsino NLSPs, there are many more final states to consider, involving 
photons, W's, Z's, and Higgses in various combinations. We will focus mostly on final states 
involving energetic photons and/or leptons, as these are clean experimental objects with 
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Figure 2: Cross sections for bino, wino and Higgsino NLSPs, as a function of the chargino mass. In the 
case of binos, the cross section shown is for production of a wino state above the bino, because direct bino 
production is negligible. Shown also are the rates for various final states, obtained by multiplying the cross 
sections by the relevant branching ratios from section [2] In order to not clutter the plots, we have not 
included the rates for mixed final states, as these arc generally small and will be irrelevant to our analysis. 

small, well-characterized backgrounds. The one exception is for Higgs-rich Higgsino NLSPs, 
where we will also consider all-hadronic final states consisting of some number of 6-jets. 

For the case of wino co-NLSPs, production of XiXi leads to final states with W ±/ y + Fj t 
and W ± Z + J£ T , while XiXi gives W + W~ + J£ T . (Direct production of XiXi is negligible, 
since this proceeds only through mixing with the heavy Higgsinos.) The following simple 
S/\fB "rules of thumb" indicate that only the first is a viable search channel for wino 
co-NLSPs: 
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W(£i')'~f+ Fj t . Here the irreducible SM background is W[£v) + 7. The cross section for 
this is a rapidly falling function of photon p?; for Pt(j) > 20 GeV, one finds a ~ 1500 
fb at the Tevatron. As can be seen from the upper right plot in fig. [2j after multiplying 
by Br(W / — > £v) m 20%, the signal rate in this final state ranges from 200-10 fb for 
ttinlsp between 100-200 GeV. So S/y/B > 1 for a sizeable range of NLSP mass, and 
this is a promising channel in which to search for wino co-NLSPs. 

W(£u)Z(££) + Fj T . Here the irreducible background is simply SM WZ production, 
whose rate in the trilepton plus J£ T final state is ~ 30 fb at the Tevatron. Meanwhile, 
fig. [2] shows that the signal rate in the W(£u)Z(££) + J£ T final state peaks at ~ 1 fb 
around ttlmlsp ~ 140 GeV. So S/yfB is never 0(1) here, and this is not a promising 
channel. 

W (jets) Z (££) + Fj T . This final state has a signal rate that peaks at ~ 4 fb also at 
m NLSP ~ 140 GeV. There are no irreducible SM backgrounds here; the dominant 
backgrounds turn out to be Z+jets and ti. However if we factor in a signal acceptance 
that is at most 30% (from a simulation of the rj distribution and an estimate of the 
geometric acceptance), there simply will not be enough events in this final state to 
offer any sensitivity. 

W[£v)W[£v) + Fj t . The irreducible background for this final state is SM WW pro- 
duction. The rate for this in the ££ + J£ T final state is ~ 300 fb. According to fig. 
[2j the signal goes from 25-2 fb as m^LSP goes from 100-200 GeV. Factoring in a 30% 
geometric acceptance, we can barely get S/\/B ~ 2 in 10 fb" 1 . So this is also not a 
promising channel to search for wino NLSPs. 



In the case of Higgsinos, we have the three different cases discussed in section 2A photon- 
rich, Z-rich and Higgs-rich Higgsino NLSPs. The relevant final states were considered in 
[30H33] , and here we will review them all. 

For the photon-rich case (low-mass Higgsino NLSPs with moderate Mi and M2), the best 
channel is again obviously 77 + J£ T . 

For the Z-rich case, the leptonic decay of the Z boson provides us with clean final states 
with tractable backgrounds, an advantage that outweighs its small branching fraction. The 
possible final states are Z(£ + £- )Z(jets) + ^ T and Z(£ + £~)Z (£'+£'-) + T . Rough S/^B 
estimates indicate that both channels could be useful: 
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• Z (£ + £~ ') Z (jets) + Here there are essentially no irreducible SM backgrounds, only 
reducible ones such as Z + jets, ti and dibosons. The signal rate is sufficiently high, 
ranging from 30-5 fb for m^LSP between 100-200 GeV, according to the lower left plot 
in fig. [2j So this channel is a promising one. 

• Z(£ + £~)Z(£' + £'~) + Fj t . This has no SM irreducible background, but the signal rate 
is small due to the By(Z — > £ + £r) 2 suppression. It ranges from 1.5-0.2 fb for ttinlsp 
between 100-200 GeV, according to fig. [2} So the channel is clean, but whether one 
can discover anything here depends on how high the signal acceptance can be made. 

For the Higgs-rich Higgsino NLSP in the hh + Fj T final state, the lower right plot in 
fig. [2] indicates that the total cross section peaks at ~ 150 fb at m^isp ~ 140 GeV, and 
it decreases to ~ 50 fb at m NLS p ~ 200 GeV. By comparison, the SM (W/Z) + h rate is 
only ~ 70 fb at the Tevatron for ~ 115 GeV, after factoring in the leptonic branching 
fractions of the W and Z. Since this is the preferred search mode for a light Higgs at the 
Tevatron, we see that the Higgs-rich Higgsino NLSP can have comparable rates to the SM 
Higgs at the Tevatron. This offers an extremely interesting possibility to discover not only 
BSM physics, but also the light Higgs at the Tevatron! 

Here the backgrounds are much more difficult to estimate, so we will not attempt any 
S/ \/~B comparisons. Instead we will only enumerate the possible final states and estimate the 
signal rates. The final states are dictated by the decay modes of the Higgs. Assuming a light 
Higgs (m h ~ 115 GeV), the largest are [50J: Br(6 66) ~ 70%, Br(/i -> WW*) ~ 10%, 
Br(/i — > t + t~) ~ 7% (we are ignoring h — > gg). From this, it is clear that to get an 
appreciable signal rate, we can only look at inclusive final states involving 6-jets. If we 
require at least two tagged 6-jets in the final state, then using a ~ 50% 6-tagging efficiency, 
the signal rate becomes 10-30 fb. The ti cross section is 7.5 pb, and cross sections for Wbb, 
Zbb, and diboson production are also of this order. Multi-jet backgrounds with fake Fj t or 
J£ T arising from semileptonic 6 or c decays are also important. Clearly, shape information and 
not just counts will be important for reducing backgrounds. Requiring more 6-jets reduces 
the signal rate accordingly, but presumably reduces the SM background by an even larger 
factor. So perhaps this is a promising way to search for Higgs-rich Higgsino NLSPs. 

Finally, let us briefly mention the possibility of mixed final states: 7 + Z + J£ T , 7 + 6+ Fj t , 
or Z + h + Fj t . Because Br(Z — > £ + £~) rs 6% and photon channels are very clean, the 7 + Z 
channel is never preferable to 7 + 7 + J£ T in the region where the photon branching fraction 
is significant. Furthermore, it is clear that 7 + 6 can only be useful in a very small region of 
parameter space where 77 = — 1 and tan/? rs 1 (see the plots in Figure [T]). We will comment 
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briefly on this scenario in Section [7j it suffers from large QCD backgrounds. Finally, an 
inclusive search strategy for Z{— > £ + £r) + J£ T + X will do just as well for the mixed Z + h 
case as for Z(£ + £~)Z (jets) + Fj t (one could inspect additional jets for 6-tags and mass peaks, 
though in practice the reach is too limited for this to be of much interest). 

To summarize, we have seen that the best channels to search for the various types of 
neutralino NLSPs are: 77 + J£ T for bino NLSPs and photon-rich Higgsino NLSPs; W[£v) + 
7 + F, T for wino co-NLSPs; Z(££) + F, T + X and Z(££) + Z(£'£') + T for Z-rich Higgsino 
NLSPs; and multi-6-jets plus J£ T for /i-rich Higgsino NLSPs. In the following sections, we 
will apply existing Tevatron analyses in these final states to constrain general neutralino 
NLSPs and describe how these analyses can be improved to optimize the sensitivity. These 
sections are organized around the different final states, so the reader who is interested in 
only specific final states can skip to the relevant (sub)sections. 

4 Searches for 77 + p T and their Applications 

The classic search for gauge mediation is the diphoton plus missing energy signature. It is 
a clean channel, with SM backgrounds arising primarily from QCD processes (77, jj, and 
H) with fake J£ T , and from electroweak processes (W/Z + 77; and W/Z + 7 + j, W/Z + j) 
with real Fj t . This channel is primarily useful for constraining bino-like NLSPs, in which 
the branching fraction to photons is guaranteed to be large. The signal rate is essentially 
that shown in fig. |2j and the SM irreducible backgrounds are virtually nonexistent. So this 
is an extremely promising channel, where we are only limited by the signal acceptance. 

There are numerous Tevatron searches in this channel [TTl - fl5] . the most recent of which 
comes from CDF with 2.6 fb -1 integrated luminosity [13] . This search required two central 
13 GeV photons which are not back-to-back, of 200 GeV, and Fj t of 3-sigma significance. 
Zero events were found with an expected background of 1.2 events. This implies the following 
95% confidence limit on the cross section times branching ratio for general neutralino NLSPs: 

o tot x Br( 7 ) 2 x e < 1.2 fb (4.1) 

where e is the signal acceptance of the CDF analysis. It is also trivial to compute the 
projected 95% confidence limit that will be set by the CDF 77 + J£ T analysis with 10 fb -1 
of integrated luminosity: 

a tat x Br( 7 ) 2 x e < 0.55 fb (10 fb" 1 ) (4.2) 
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Figu re 3: Various examples of 95% exclusion contours (shaded regions) in the (/j,, Mi) plane, derived from 
(4.1), for the most recent 2.6 fb _1 CDF search for anomalous 77 + J5 T production. Also shown are the 
projected exclusion contours with 10 fb _1 of data (dashed lines). On the left: exclusion contours for tan/3 = 
20, rj = +1 and various values of Mi. On the right: the same, but for tan/3 = 1.5, 77 = — 1. 



One can check that with 10 fb^ 1 of data, the 77+ J£ T channel does not offer any opportunity 
for 5<T discovery. 

In fig.[3j examples of exclusion contours are shown in the (/z, Mi) plane, for various choices 
of M2, tan/3 and 77. These plots, unlike the remainder of the plots in this paper, are based 
solely on analytic calculations of the cross sections, rather than a Pythia and PGS simulation, 
because here we are able to rely on the experiment's quoted value of eS 3 ' 1 These plots 
illustrate the power of 77+ J£ T in constraining the general neutralino NLSP parameter space, 
and not just the specific simplifying limits described in the previous sections. Taking M2 = 
1000 GeV illustrates how mixed bino-Higgsinos can be excluded with 77+ J5 T . Similarly, the 
curves with M 2 = 400 show how photon-rich Higgsinos (large M 1; small fi) can be excluded. 

( 3 )To make the exclusion plots, we have assumed a constant overall acceptance of e = 7.8%, as quoted in 
the CDF analysis for the MGM benchmark model SPS8 [51 . In reality, the acceptance varies somewhat as 
we move about in the parameter space, due mostly to the dependence of the photon isolation requirements 
on the splittings between different -ino states. In our PGS simulations, we have found variations by as much 
as a factor of 2 from the efficiency quoted in the CDF study. However, because a x Br is varying rapidly 
across the parameter space, we find that even if our estimate of the efficiencies is off by a factor of 2, our 
reach estimates vary by less than 5%. 
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Finally, with M 2 = 250, we see that 77 + Fj t can also exclude mixed wino-Higgsino NLSPs 
(e.g. M2 = 250 and \i w 200 with M\ — ► 00), where the branching fraction to photons is 
due to partly to the wino content, as in |2.7 In the last two examples, note the difference 



between the left and right plots in fig. [3] - in the latter, the bounds are weaker because the 
branching fraction to photons via the Higgsino is tan j3 suppressed. 

For wino co-NLSPs and heavier Higgsino NLSPs, the 77 + F^ T channel is not useful. For 
Higgsino NLSPs, the branching fraction to photons is a result of phase-space suppression in 
the other channels, so making the NLSP heavier turns off the rate in 77 + Fj t . For wino 
co-NLSPs, the 77 + J£ T signature can only come from direct production of XiXx- However, 
as noted above, the cross section for this is negligible in the wino co-NLSP limit, since only 
Higgsino neutralinos can be directly produced. 

Finally, let us specialize to the case of bino NLSPs and discuss the bound from 77 + Fj T 
here. As discussed in the previous section, in the case of a bino NLSP, the production 
cross section is determined by the mass of the heavier chargino states and is essentially 
independent of the bino mass. Here we consider the case that these heavier states are winos 
(i.e. \M 2 \ <C Comparing (4.1) and the plot of the -ino cross section in fig. |2j we see 



that the lower bound on the chargino mass (M 2 ) is( 4 ) 

m x ± > 270 GeV (4.3) 

In 10 fb _1 , the projected bound on the chargino mass is: 

m x ± > 300 GeV (10 ftT 1 projected) (4.4) 

We should point out that the kinematics and event topologies can be very different across 
parameter space. So if an excess is observed in this channel in future searches, such differences 
can be important in characterizing the physics. The main point is that for bino NLSPs, the 
rate is set by the masses of heavier winos, which can be far above the bino mass. Thus 
events will typically have energetic decay products coming from W — > B + X . For instance, 
in MGM, the decay chain can go through a slepton and one expects energetic leptons or 
taus in the events. However, for Higgsino NLSPs, the production cross section is determined 
mostly by the mass of the NLSP itself, with charged states split from neutral states by an 
amount of order 1-10 GeV. The decay products of heavier -inos down to the NLSP will be 



( 4 )This result is valid in the large /i limit, and is the weakest possible bound in the scenario /i 3> M2 > Mi. 
The stronger bound of 300 GeV quoted by CDF relics on a particular SPS benchmark model point. 
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Figure 4: Feynman diagrams for two-photon processes arising from neutralino and chargino production. At 
left, the typical process in MGM, where xt are mostly wino and decay through sleptons to the mostly-bino 
X\- The final state includes energetic tau leptons. At right, a typical process with mostly-Higgsino NLSPs, 
which are produced directly. The small splitting between xt an d Xi leads to a three-body decay through 
off-shell W with very little phase space, so there are relatively soft leptons or jets in the final state. 

much softer, and the event could contain little additional activity beyond 77+ E T . Examples 
of the different decay chains are shown in fig. |4} 

Various other Tevatron searches involving energetic photons and missing exist [52T 
154"] . We have analyzed them in some detail; while some parts of parameter space can be 
excluded with these results, we find that the limit from 77 + E T is always much stronger, 
and so we will not discuss them in detail. 

5 Searches Relevant to wino co-NLSPs 

5.1 Searches for 7 + W + J£ T 

CDF has published a search for 7 + t + E T with 0.93 fb -1 of data [55]. They selected for 
at least one isolated photon and at least one isolated lepton (e or \i) with pt > 25 GeV 
and \t]\ < 1. They also required E T > 25 GeV. They found 163 events with an expected 
background of 150.6 ± 13.0. This null result sets a 95% confidence limit on the cross section 
times branching fraction for general neutralino NLSPs: 

a x Br x e < 40 fb (5.1) 

With 10 fb -1 , the projected bound is 

a x Br x e < 8 fb (5.2) 
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Figure 5: Plots of a x Br x e (i.e. the number of events per fb _1 ) in the £ + 7 + Fj T channel coming from wino 
co-NLSPs. The left plot shows the number with J£ T > 25 GeV, corresponding to the cut used in the CDF 
study. The right plot shows the number with Fj T > 50 GeV. The solid lines indicate the 95% confidence 
limit calculated from the CDF study, while the dashed lines indicate the projected bound in 10 fb _1 . 

These bounds can be used to constrain wino co-NLSPs, in which one side is a chargino that 
decays promptly to W + G, and the other side is a neutralino that decays to 7 + G. 

Shown in fig. [5] is a plot of a x Br x e in the 7 + I + J£ T channel coming from wino 
co-NLSPs. Here the signal acceptance e was computed from PGS simulations, as described 



in section 1.2 From the figure, we see that the existing CDF search already bounds wino 
co-NLSPs, 

m NLSP > 115 GeV (5.3) 

With 10 fb _1 , the bound potentially becomes 

m NLS P > 155 GeV (10 fb^ 1 projected) (5.4) 

With simple modifications to the existing analysis, we find that the sensitivity can be 
extended significantly. In particular, we find that the Fj T coming from SM WW production 
is generally much less than the Fj t in the signal. The second plot in fig. [5] shows our estimate 
of the bounds if a Fj t > 50 GeV cut is applied: 

m NLSP > 135 GeV (£ T > 50 GeV) (5.5) 
m NLSP > 170 GeV (£ T > 50 GeV, 10 fb" 1 projected) (5.6) 
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The backgrounds were derived by using the J£ T distribution shown in fig. 4 of [55]. We see 
that a harder J? T cut improves the existing bound and the projected bound by 15-20 GeV. It 
is possible that an even harder J£ T cut would improve things further, but for that we would 
need a more accurate determination of the backgrounds at very large Fj t than is available 
in [55]. Also, it would be interesting to investigate cutting on other kinematic variables such 
as H T and M T .& 

Given the significant statistical improvement that a harder Fj t cut and 10 x more data 
will bring to the existing search, it is also worth mentioning the discovery potential in this 
channel. We estimate that a 5a discovery is possible up to rriNLSP ~ 140 GeV, and 3a 
"evidence" is possible up to m NLS p ~ 160 GeV. 

Finally, let us briefly discuss the possibility of using 7 + jets + J£ T to constrain the wino 
NLSP scenario. Here obviously one is substantially increasing the signal rate through the 
higher branching fraction of W — > jets, but one is also increasing the SM backgrounds by 
introducing hadronic activity into the event. There do not appear to be any public Tevatron 
searches in this channel; the closest seems to be [52], which unfortunately searched specifically 
for events with displaced photons. Also, preliminary results on 7+jets(+Fj T ) were presented 
in [56]; however, at present this analysis has no Fj t cut, and therefore the backgrounds are 
simply too high for it to be useful. Without a more detailed analysis, it remains to be seen 
whether signal or background will win out in the end. 

5.2 Other searches 

There have been several Tevatron analyses of the W + (£ + u) + W~{l~v) final state, including 
a CDF measurement of the W + W~ cross section with 3.6 fb _1 of data |57| . and a D0 mea- 
surement of the same cross section with 1 fb _1 of data [58J. We find that neither search 
constrains wino NLSPs even with 10 fb _1 of integrated luminosity. So this channel does 
much worse than the 7 + W{tv) channel considered in the previous subsection. The reasons 
appear to be: (1) the relative signal suppression due to the extra factor of Br(W — > hi) and 
the fact that cr(pp — > XiXi) ~ \ a ipP ~^ XiXi) i n the range of NLSP masses of interest; 
and (2) the high background rate ~ O(10 2 ) fb from a variety of sources. In particular, as 
discussed in section |3j SM WW production has a rather high cross section relative to the 
signal, and it is very difficult to eliminate this source of background. We have investigated a 
number of simple kinematic variables, such as Fj T , and have not found one that can signifi- 

( 5 >The Standard Model W-y process has a large charge asymmetry, which may also be useful in distinguish- 
ing signal from background in this channel. We thank Yuri Gershtcin for this observation. 
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cantly reduce SM WW while keeping enough signal. Perhaps a more sophisticated kinematic 
quantity such as Mt,2 can help [59j ED] ■ 

We have similarly examined the CDF analysis in the Z(e + e~) + W^jets) + J£ T final state 
using 2.68 fb -1 of data [61]. Here the main analysis cuts were: central, high-pr, opposite- 
sign di-electron pair forming a Z; a dijet pair forming a W; and J£ T >40, 50 or 60 GeV. 
With these three different J£ T cuts, CDF observed 7, 2, and 1 event with 6.41, 3.76 and 2.02 
expected, respectively. Simulating the signal, we find that this analysis does not constrain 
wino NLSPs even with 10 fb _1 of data. The problem again seems to be the low signal rate; 
as discussed in section |3j the wino NLSP cross section is below 4 fb everywhere, and the 
signal acceptance is at most 30% just from geometric considerations. 

6 Searches Relevant to Z-rich Higgsino NLSPs 

6.1 Searches for Z {£+£-) + ^ T + X 

Results in this channel include a CDF search for Z{£ + £~) + J£ T (+X) [62], a D0 measurement 
of the SM Z(£ + t) + Z(vv) cross section [63], and a CDF search for Z(e + e~) + VK(jets) + Fj t 
|61j . The dominant SM backgrounds for all these Tevatron searches consist of Z + jets, 
W + jets, ti, and dibosons. 

We find that none of the existing studies currently rule out any part of the general 
parameter space. However, with 10 fb^ 1 of data, all of these searches have the potential to 
exclude regions of Higgsino and wino-like parameter space, and some, with cuts optimized 
for our scenario, even have the possibility for genuine 5<r discovery. In particular, the search 
for Z(£ + £~) + Fjrp + Hj* in [62] appears to be the most effective. So let us now briefly describe 
the search and its implications for general neutralino NLSPs. 

The requirements were pairs of central, opposite-sign electrons or muons (one with > 
20 GeV, one with p T > 12 GeV) reconstructing a Z within the window 66 GeV < M inv < 
116 GeV. Further requirements were Fj t > 25 GeV with 3-sigma significance (i.e. MetSig 
> 3, as defined in [62]) and Ht > 300 GeV. With these cuts, CDF observed 6 e + e~ events 
with 6.6 ± 1.0 expected, and 6 events with 3.5 ± 0.5 expected, with an integrated 

luminosity of 0.94 fb -1 . 

The projected 10 fb -1 exclusion contours in the (/i, Mi) plane are summarized by the 
plots in fig. [6] For comparison, the 77 + Fj t contours for the same choices of the parameters 
are also shown. We see that the two channels are mostly complementary in their coverage 
of the parameter space. The exclusion contours from the Z + Ht + J£t search are largely 
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Figure 6: Various examples of projected 95% exclusion contours (blue, dashed) in the (/1, Mi) plane, based 
on the CDF search for anomalous Z(£ + £~) + Ht + Fj t production. The contours are derived assuming 10 
fb _1 of data. Also shown are the projected exclusion contours (green, dotted) for the 77 + J? T search, for 
the same choices of the parameters. 

independent of the other parameters (M2, tan (3, 77), and the typical situation is characterized 
by the left plot. The one exception, shown in the right plot, is when tan (5 is small and rj = — 1. 



As discussed in section 2.4 in this case the branching fraction to Higgses takes over at larger 
NLSP mass, and this weakens the Z + Ht + sensitivity. 

The pure Higgsino NLSP limit corresponds to Mi — > 00 in the plots in fig. [6j We see 
that in this limit, the projected bound is 

m NLSP > 150 GeV (10 ftT 1 projected) (6.1) 

essentially independent of the other parameters, except again when 77 = — 1 and tan/3 — > 1. 
In this case, we estimate that the Z(£ + £~) + Ht + search should be able to cover the 
entire Z-rich range of NLSP mass. For instance, if rj = —1 and tan/3 = 1.5 (the case shown 
in figs. [l]and[6]), we find that the projected bound is tunlsp > 130 GeV. For larger values 
of mjsiLSP the branching fraction to Higgses takes over. 



23 



0.5 - 




100 200 300 400 500 600 50 100 150 200 250 



HT (GeV) MET (GeV) 



Figure 7: Comparison of distributions for a Higgsino NLSP signal point (// = 170 GeV, Mi, M2 — > 00, tan (3 
= 20) and backgrounds. Signal and total background are independently normalized to unit area. In every 
case we plot only for those events passing the cuts in the Z(e + e~) + J5 T analysis of [62 , i.e. J? T > 25 GeV, 
MetSig> 3, and a central Z decaying to e + e~. For the Mi nv (jj) distribution, the invariant mass of the 
hardest two jets is plotted, and the zero bin is populated by events with less than two jets. 

6.2 Suggestions for Optimizing Z + J£ x + X Searches 

We have seen in the previous subsection that with enough integrated luminosity, existing 
searches at the Tevatron can potentially exclude Higgsino NLSPs with mass up to ~ 150 
GeV. However, these searches are not targeted specifically to the Higgsino NLSP scenario, 
and as we will discuss in this subsection, some simple optimizing of the cuts can easily extend 
the reach by ~ 20 GeV. 

In fig. [7j we compare a few kinematic distributions for the main SM backgrounds and 
a particular Higgsino NLSP signal point. These distributions were chosen to highlight the 
differences between background and signal. And apart from M inv (jj), they also happen 
to be the kinematic variables used in the CDF analysis. While we have not performed a 
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comprehensive analysis, other distributions that we have examined, such as Nj ets , Nb-j e ts, 
cos A(fi(Z, J£ T ), etc., do not seem to separate background and signal as cleanly. 

From the figure, we see that simple modifications of the CDF cuts can potentially improve 
the sensitivity to our signal hypothesis. For instance, tightening the Z mass window cut may 
be useful. This was set at 66 GeV < Mg+i- < 116 GeV in the CDF analysis [62], so as to 
include the entire range in which the Z peak in Drell-Yan rises above background. We see 
from the figure that a narrower range, such as M^+£- e (80, 100) GeV, will still keep most of 
our signal, while cutting down on backgrounds without real Zs (W + W~ and tt) significantly. 

The Ht distribution in fig. [7] shows that the Ht > 300 GeV cut used in the CDF analysis 
tends to eliminate too much of the signal. This is simply because the Ht in the signal is 
coming entirely from the pair-produced -ino decay products, and the -ino masses considered 
here are in the ~ 100 — 200 GeV range. According to the figure, a lower Ht threshold, such 
as Ht > 200 GeV, would boost the signal acceptance while still eliminating most of the 
backgrounds except tt. 

Finally, the Fj T distribution in fig. [7] shows that a harder J£ T cut, such as Fj T > 40 GeV, 
can improve S/ VB significantly. 

Motivated by these considerations, we find by optimizing the analysis of Z + J£ T + X 
along these lines, the reach can be improved to 

m NLS P > 170 GeV (10 fb -1 projected) (6.2) 

for the Z-rich Higgsino NLSP case. Also, we find that 5a discovery is possible up to rriNLSP ~ 
135 GeV and 3a evidence up to m NLS p ~ 160 GeV. The background and signal (ttinlsp 
1 70. 7] = +1 and tan/3 = 20) counts are shown in table [I] for the following set of cuts: 80 
GeV < M l+l - < 100 GeV, H T > 200 GeV, T > 40 GeV, MetSig> 3. 

Additional cuts could potentially play a role, but we have not investigated them in detail. 
For instance, tt background can be reduced by a cut on the second-hardest jet in the event. 
Mass window cuts could be used to attempt to isolate a second Z boson decaying to jets (as 
indicated in the M inv (jj) distribution shown in fig. J7j), or a Higgs decaying to bb (perhaps 
with 6-tagging, as well). As we discussed in the previous section, the splitting between 

( 6 ^A brief comment on our background simulations: with a harder J5 T cut, we are potentially underesti- 
mating Z+jets, which comprises ~ 20% of the background after the cuts described here. The reason is that 
PGS does not accurately reproduce the tails of the Z + jets JC T distribution, since it is mostly coming from 
mismeasurement. However, our use of the MetSig variable defined in [62 a should help, since this involves 
detailed information about jet energy resolution. Indeed, our simulations are able to reproduce to 10-20% 
accuracy all of the analyses in [55] which involve MetSig. So hopefully this means that the background 
estimates used in the modest extrapolation described here are not too far off. 
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Table 1 : Rates for a proposed Z + + J£t search 



N b9 (e + e-) (fb) 


N bg (fi + fx-) (fb) 


N sig (e + e-) (fb) 


iV si9 (/i+/i-) (fb) 


7.4 


5.8 


1.4 


1.1 



charged and neutral Higgsinos can be quite small, and will lead to fairly soft jets or leptons. 
So selecting events with a soft, isolated lepton from the decay xt ~~ > ^ +Z/ X? could help to 
purify the signal. Needless to say, if an excess is observed in this channel, any of these 
features, even if not useful for obtaining good results in a counting experiment, could help 
in understanding whether or not an NLSP explanation fits the data. Further improvement 
might come from considering Z — > t + t~ final states, but due to the neutrinos in the r 
decays, such events would be unlikely to pass the tight Z mass- window cut, and one would 
have to reconsider how to reject backgrounds. 

Finally, we would like to emphasize that with very simple, naive cuts alone we have seen 
that the sensitivity in this channel can be extended by ~ 20 GeV. After these cuts, we are 
still finding ~ 25 signal events in 10 fb -1 . This strongly suggests that with smarter cuts 
and more sophisticated experimental methods (e.g. neural nets), the sensitivity could be 
extended even further. 

6.3 Searches For Z(£+£~) + Z {£'+£'-) 

Here we discuss Tevatron measurements of the SM ZZ cross section using the Z(£ + £~) + 
Z{£' + £'~) final state. The most recent published analyses are one from CDF with 1.9 fb -1 
of data |64J and one from D0 with 1.7 fb -1 of data |65| . In addition, there is a very recent 
public note of a CDF analysis with 4.8 fb^ 1 of data [66J. All of these analyses require two 
opposite sign dilepton pairs satisfying various kinematic and isolation requirements. 

We have chosen to reproduce the published 1.9 fb -1 CDF analysis for simplicity. They 
observe 2 events in 1.9 fb -1 , and the SM prediction (coming almost entirely from ZZ pro- 
duction) is 2.004+°;°^ ± 0.210. 

Simulating the signal, we find that the existing CDF analysis has no sensitivity in the 
general neutralino NLSP parameter space. Furthermore, even in 10 fb -1 the projected 95% 
exclusion contours are negligible, and are completely covered by 77 + J£ T and Z + Ht + J£ T . 
The reason for this is simply because the overall rate is suppressed by two factors of the small 
branching fraction of Z — > £ + £~, and so S/\rB < 1 everywhere in the parameter space. 
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Z(/ + r)+Z(f + f ~)+(MET>25 GeV) 




^nlsp (GeV) 

Figure 8: The projected number of Higgsino NLSP events in 10 fb _1 in a hypothetical search for Z(£ + £~) + 
Z{t' + l'-) + £ T) with F, T > 25 GeV. 

However, things could be improved with a J£ T cut, since the background (SM ZZ pro- 
duction) should have very little intrinsic Fj t . We do not attempt to simulate fake Fj t from 
SM ZZ production, we simply assume that a reasonably sized J£ T cut, e.g. Fj t > 25 GeV, 
should reduce the background to JV& <C 1. Then we would like to have at least a few events 
to claim a discovery. A plot of the estimated number of Higgsino NLSP events in 10 fb _1 in 
this hypothetical search is shown in fig. |8j Here Mi, M 2 are taken to be very large, \x > 0, 
and tan/3 = 20. We see that, assuming the SM background is very small, we could conceiv- 
ably have a "discovery" for uinlsp up to ~ 140 GeV. Of course, given the low statistics, a 
true estimate of the discovery potential would require a much more accurate estimate of the 
small SM background. We should also emphasize that the low statistics highly motivates a 
DO/CDF combined analysis in this channel. 

7 Searches Relevant to /i-rich Higgsino NLSP 

Finally, we come to the most difficult scenario to analyze, in which the NLSP decays dom- 
inantly to h + G. It occurs in a relatively small part of parameter space (Higgsino NLSP 
with 7] = —1 and tan/3 ~ 1), but is interesting, especially for the prospect of observing the 
Higgs at the same time as physics beyond the Standard Model. Because one can have Higgs 
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bosons on both sides of the event, each decaying to bb, the signals can involve anywhere from 
two to four 6-tagged jets along with missing Et- 

Before considering the multi-fe channel in detail, we note that the channel 7 + b + j '• + E T 
is potentially useful in a small slice of parameter space where tan (3 « 1 and Mi « — |//. (In 
general, one does not have -Br (7) ~ Br(h) ~ 0(1); see the lower right-hand plot of Figure 
[IJ) Early in Run II, a Tevatron SUSY Working Group considered this channel in some detail 
[26] . choosing a "Higgsino Model Line I" in this slice. Recent observations in this channel 
show a rate of ~15 fb for events with E T > 50 GeV, consistent with the backgrounds, which 
are fake-dominated, with large systematic uncertainties [HZj. The early working group study 
relied on lower background estimates extrapolated from the small Run I dataset. Hence, we 
do not consider 7 + b + j + E T a very promising channel. We will confine our attention to 
studies involving multiple b jets plus E T in the remainder of this section. 

Backgrounds in the multi-6+E T channel are large and uncertain. We can get a sense of 
how feasible a search is by looking for similar existing new physics searches. Searches for 
more than 2 fe-tagged jets are motivated by the prospect of MSSM Higgs bosons [68, 69J. 
Most recently, such a study from D0 JBH] required 3 jets with tight 6-tags in events with up 
to 5 jets and loose acceptance cuts. In a sample of 1 fb -1 , they found 3224 3-jet events, 2503 
4-jet events, and 704 5-jet events. These samples are much too large for our signal to show 
up in a simple counting experiment. Furthermore, D0 set limits using only the shapes, not 
counts, due to uncertain backgrounds. The M{bb) distributions in the data peak just above 
100 GeV, so a cut in this variable will not help us. A hard E T cut is more promising, but 
we do not have access to the E T distribution in the data to read off a reliable count. 

We learn more from existing searches for bb+ E T . A CDF search for ZH/WH production 
[70] and a D0 search for leptoquarks or sbottoms [71] each look for two 6-jets and E T , but 
due to a veto on events with 4 or more jets, they fail to constrain the Higgsino scenario. 
However, similar searches are potentially constraining, as we will see shortly. More interesting 
is a CDF search for gluino-mediated sbottom production using 2.5 fb -1 [72], which looks for 
E T > 70 GeV and two 6-tagged jets with some basic kinematic cuts. In a sample that further 
rejects leptons and requires E T to be far from jets in 0, they observe 451 events with 506 
± 144 expected. Note that the rates and the uncertainties are large, even with a strict E T 
cut. Top quark and heavy-flavor multijet backgrounds are dominant. 

Owing to their high multiplicity and their all-hadronic, heavy-flavor content, we do not 
attempt to simulate the multi-6+E T backgrounds in detail, as in the previous sections. How- 
ever, we can give some preliminary remarks. We have simulated ti and diboson (W 7 W 7 , WZ, 
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and ZZ) samples in Pythia (37], and Zbb, Wbb, and bbbb samples in MadGraph [38]. These 
background samples were processed with PGS, which includes a parametrization of CDF's 
SecVtx 6-tagging efficiency as a function of Et and 77. For a messy hadronic channel like this 
one, backgrounds should be simulated with matching, and run through the collaboration's 
full detector simulation, to obtain realistic estimates. Nonetheless, we offer a rough guide to 
how to begin to design a search. 

We first examine a set of cuts that asks for 50 GeV of Fj t , at least two 6-tagged jets, 
and no central electrons or muons. The rate is still very large (as suggested by the results of 
Ref. [72]), so we further require that at least one pair of 6-tagged jets has an invariant mass 
satisfying 60 GeV < M{b X: b 2 ) < 200 GeV, with AR(b 1 ,b 2 ) < 2.5. The latter requirement is 
imposed because 6-jet pairs in the background samples are unlikely to have large invariant 
mass if they are not back-to-back. (See the right-hand plot in Figure [9j which we will discuss 
in more detail shortly.) The results are tabulated in Column A of Table [2j 

We next consider a second set of events with more fe-tags: either three or more tight b- 
tags, or two tight plus at least one loose 6-tags. This substantially reduces the ti background, 
which rarely has three 6-tags unless one of the W decays gives a charm quark (which has 
~ 10% probability of being 6-tagged). Demanding more 6-tags allows us to relax the J£ T cut 
somewhat, so we only ask for 40 GeV of Fj t in this sample. The results are in Column B of 
Table U Note that now the bbbb sample is much more important than it was in the sample 
with two tags and a harder Fj t cut. 

The final column in the table displays the rate for events with four 6-tagged jets (up to 
two of which can have only loose tags), all with pt > 15 GeV and Fj T > 30 GeV. Unlike in 
columns A and B, no cut on the invariant mass of pairs has been imposed. The background 
rates, especially for bbbb, rise steeply when the pr cut on all jets, or the J£ T cut, is loosened. 

From the table, one can see that if these samples encompass all the important back- 
grounds, and if the 6-tagging efficiencies are modeled well by PGS, and if the systematic 
errors are small enough to be neglected, then a counting experiment would be able to ex- 
clude some region from fi « —150 GeV to \i ~ —180 GeV at 95% confidence. (We note 
that below fi ~ 150 GeV, the Higgsino NLSP mostly decays to Zs, and this is covered by 
the analyses described in the previous section.) In practice, we expect that this is not the 
case. In particular, our samples overestimate the ti numbers reported in Ref. [72] and un- 
derestimate the QCD heavy-flavor contribution (because we have not included bbjj, etc.), 
and the estimates reported there carried large systematic uncertainties, which we have not 
even attempted to quantify for our samples. It is clear that a search will have to rely on 
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Table 2: Multi-6 + ^ x rates in fb 



Sample 


A. bb + F, T >50 


B. bbb + F, T >40 


C. bbbb+fi T >30 


it 


77.2 


16.8 


1.7 


Wbb 


12.4 


1.4 


0.0 


Zbb 


6.1 


0.8 


0.1 


bbbb 


4.8 


9.1 


1.6 


Diboson 


2.7 


0.2 


0.0 


Total a 


103.2 


28.3 


3.4 


m NLSP = 140 GeV 6 


5.2 


3.1 


1.0 


ttinlsp = 160 GeV 


7.1 


4.1 


0.9 


m NLSP = 180 GeV 


6.2 


3.3 


0.7 


m NLSP = 200 GeV 


4.5 


2.3 


0.5 



a Includes only simulated backgrounds - not comprehensive. 

6 In all signal points listed, Mi = 500 GeV, M 2 = 800 GeV, 77 = — 1 and tan/3 = 1.5. 



the shapes of distributions, not just counting. From the results in Refs. [70tf72] . we can 
be fairly confident that neural nets or other statistical techniques can eliminate much of 
the QCD multijet backgrounds with fake J£ T or Fj t arising from semileptonic b or c decays 
(which will tend to be aligned with the jet direction), so the results in the table may not 
be wildly optimistic. It appears that the largest background, especially if only two 6-tags 
are demanded, will be tt. Statistical analysis might be able to label events as "top-like" and 
improve the discrimination. In the end, any limit set in this channel will probably rely on 
the fact that tt events have a broad distribution of M(bb), which is more likely to be near 
rrih in the signal. 

One might ask whether the signal can be isolated just from invariant mass distributions. 
For instance, one could apply only a loose J£ T cut, demand four 6-tagged jets, and look for a 
pairing where M(bib 2 ) ~ M(6 3 &4), with signal accumulating near the Higgs mass. However, 
we find that in the PGS output, such pairs do not cluster strongly near (m h ,m h ), and offer 
little hope of finding a peak above the background (see the left-hand plot in Figure [9]). 
Column C of Table [2] demonstrates that any technique that relies on tagging four 6-jets will 
have limited signal rate, though with a Fj t cut a signal-to-background ratio of order 1 may 
be achievable. One would want to rely on shape information, such as that illustrated in the 
right-hand plot of Figure [9j to gain confidence that a signal is real. This figure illustrates 
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Figure 9: At left: nearest pairs of M(bb) in partitions of 4 b-tagged jets, in 100 signal events (blue) and 100 
bbbb background events (red). At right: all combinatoric possibilities for invariant mass and AR between 
pairs of b-tagged jets in 50 signal events with 4 b-tags (blue) and 50 bbbb background events with 4 b-tags 
(red). The signal point is fi = -140 GeV, Mi = 510 GeV, tan/3 = 1.5, M 2 = 800 GeV. 

the correlation that was used in choosing cuts for columns A and B of Table [2] In the 
background events, pairs of jets with large invariant mass are rare, and typically such paired 
jets are nearly back-to-back (AR ^ 7r). On the other hand, for the signal events, the two 
b jets from a Higgs decay are typically not back-to-back (though they are well-separated). 
Notice that for 1 ^ AR ^ 2, M(bb) in the signal cuts off at around m^. Thus this shape 
offers some hope of discriminating signal from background and obtaining an estimate of m^. 

In practice, we expect that in an all-hadronic channel the understanding of the back- 
ground shapes and normalization will have to be taken from data rather than Monte Carlo, 
and more sophisticated techniques such as neural nets or boosted decision trees could be 
useful. Further improvement could come from D0, where 6-tagging is efficient at larger val- 
ues of r] than at CDF. (At CDF, soft lepton tagging in addition to SecVtx might also yield 
improvements beyond the efficiencies parametrized by PGS.) We expect that the general 
picture we have sketched — using a J£ x cut to reduce the overall size of the background, 
then studying the shape in the (M(bb), AR) plane to identify the signal contribution - 
could survive these complications and prove to be useful. 

Our results are broadly consistent with those of the SUSY Working Group operating 
at the start of Run II [26J. Their assumption (given the more limited Monte Carlo tools 
at the time) was that backgrounds were dominated by ti. This is reasonable, so long as 
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one only demands two 6-tags. We have found that the multijet contributions that were not 
included in early studies are important when three or more 6-tags are included, however. A 
full experimental study will, of course, also have to include bbcc, cccc, bbjj, etc. 

While such a search is challenging, and is necessary to find NLSPs only in a limited region 
of parameter space, Br(x5 - ► h + G) is generically order-one whenever m^o is sufficiently 
far above m^. Thus, this search is worth pursuing for the prospect of finding the Higgs, 
especially if other channels like Z(£ + £~) + Fj t begin to show an excess. Finding decays 
to Higgs would also be an important step in establishing the Higgsino-like nature of any 
prospective NLSP candidate. 

8 Conclusions 

8.1 Summary of Results 

In this paper, we have analyzed in detail the signatures of promptly-decaying, general neu- 
tralino NLSPs at the Tevatron. Through both simple S/ \fB considerations and detailed 
simulations, we have derived the most promising Tevatron search channels for each type of 
general neutralino NLSP: 77 + F^ T , for bino or very light Higgsino NLSP; W(£v) + 7 + J£ T , 
for wino co-NLSPs; Z(ll) + F, T + X and possibly Z{U) + Z(£'£') + F, T for Z-rich Higgsinos; 
and possibly multi-&+ Fj t for Higgs- rich Higgsinos. 

The estimated limits from existing Tevatron searches, along with the projected limits in 
10 fb" 1 for optimized searches, are presented in Table |3]for the cases of pure bino, wino, and 
Higgsino NLSPs. More generally, there are limits on mixed NLSPs arising from 77 + Fj t 
which cannot be easily summarized in the table, but are displayed in Figure |3j Note that 
for bino NLSPs, the result is expressed as a limit on the lightest chargino mass, while in the 
other cases it is a limit on m NLS p- However, in these other cases, m NLSP f& m~±, so we can 
always think of the limit as being on the lightest chargino mass. 

Somewhat surprisingly, we also saw that the Tevatron still has an opportunity for 5cr 
discovery of wino and Higgsino NLSPs. We found that in the optimized W(£u) + 7+ J£ T and 
Z(£ + £~) + Hj< + Fjrp channels (and perhaps in the very clean Z{-+ £ + £~)Z(-^ £'+£'-) + F, T 
channel), discovery is possible for wino and Higgsino NLSPs, respectively, up to mNLSP ~ 
135 - 140 GeV. 

Mostly, the role of the Tevatron will be to exclude regions of parameter space, or more 
optimistically to build up a sample of interesting statistically significant events that falls 
short of discovery. Because an NLSP that is a good target for one search (e.g. Z + J£ T + X) 
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Table 3: Summary of Results. All numbers are in GeV, and unless otherwise specified, refer to the 
NLSP mass. 



NLSP Scenario 


Search Channel 


Curr. Lim. 


Proj. Lim. 


Proj. 3a 


bino 


77+ 


?V > 270 


m~± >300 

Xi 


None 


wino co-NLSP 


W{iv) + 7 + 


135 


170 


160 


Z-rich Higgsino 


z(e+e-) + p T + x 


None 


170 


160 




Z{£ + £~)Z {£'+£'-) + F, T 


None 


160 (?) 


140 (?) 


Higgs-rich Higgsino 


multi-6+p T 


None 


160 (?) 


(?) 



has small branching fractions to other decay modes, the Tevatron will have limited ability 
to correlate interesting excesses in more than one channel. On the other hand, the events 
in a given channel may contain hints about the underlying physics that can be exploited in 
future searches at the LHC. For instance, as we have discussed, a 77 + Fj t excess with only 
soft jets in the rest of the event is suggestive of a light Higgsino NLSP, which will in general 
have a nonnegligible branching fraction to Z + G. At the Tevatron we find that 77 + Fj t 
and Z + + X have little overlap in their reach, but perhaps the LHC can fare better in 
substantiating the Higgsino interpretation with searches involving Z bosons, especially if it 
can produce SUSY events with a strong cross section. 

8.2 Future Directions 

In this paper we have focused on the specific scenario of prompt neutralino NLSPs at the 
Tevatron. As discussed in the introduction, searching for this scenario is one of the most 
well-motivated ways that the Tevatron can look for gauge mediation. However, there are 
other studies related to gauge mediation at the Tevatron, that could be interesting as well. 
We have not exploited the possibility that additional supersymmetric states, such as sleptons 
or colored sparticles, could increase the reach for gauge mediation at the Tevatron. The phe- 
nomenology of including these states in combination with a general neutralino NLSP might 
also reveal new blind spots in the Tevatron's search for supersymmetry. For instance, a Higgs- 
rich Higgsino NLSP produced in the decay of a heavier colored sparticle may avoid standard 
constraints from jets + Fj t searches. There are also extensions of the general neutralino 
NLSP scenario that do not rely on additional sparticle states. The most straightforward 
extension is to consider delayed neutralino NLSP decays rather than prompt. In gauge me- 
diation models this corresponds simply to raising the fundamental scale of supersymmetry 
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breaking without having to change anything else. Using the tracking capabilities at both 
Tevatron experiments, and the EM timing capabilities at CDF or EM pointing capabilities 
at D0, the Tevatron can set interesting limits on delayed general neutralino NLSPs [3"5] . 

While the focus of this paper has been the Tevatron, we don't intend to downplay the 
LHC's capabilities. It is natural to ask how well the LHC can do, and how quickly it will 
surpass the capabilities of the Tevatron. Given that the LHC has not yet collided beams, it 
may be quite some time before it will have delivered comparable amounts of data. 

Existing LHC studies are largely focused on the minimal gauge mediation scenario 
(though see Ref. [IH])- For instance, a recent ATLAS study [IB] of 77 + E T shows that 
this can be a very good search channel with even small amounts of luminosity. Recently AT- 
LAS has considered the case of Higgsino NLSPs in the Z(— > £ + £~) , y channel |73j . estimating 
a reach of about 135 GeV with 3 fb" 1 at -/a = 10 TeV. We have seen that the Tevatron can 
already exclude this region of parameter space. 

The main advantage of the LHC over the Tevatron would be the ability to make gluinos 
and squarks; if strongly-produced particles are accessible, there is a capacity for a much 
higher rate. However, we should not be too wedded to particular ideas about the SUSY 
spectrum. The gluinos and squarks could be sufficiently heavy that this is not a strong 
advantage' 7 '. For some channels, like searches involving h — > bb, QCD backgrounds at the 
LHC can be more severe; the Tevatron might be a slightly cleaner environment (though, 
naturally, any channel in which QCD backgrounds are dominant will be an uphill struggle 
to achieve statistically significant results at any hadron collider). 

The issues of backgrounds at the LHC and a detailed comparison of LHC and Tevatron 
reach deserve further study. We plan to examine these and other issues in detail in a 
forthcoming paper |34j . In particular, gauge mediation signatures at the LHC along the lines 
of the Z + E T + X and multi-6+ E T signatures that we discussed in this paper have received 
little attention so far. It would be interesting to consider not only simple extrapolation of the 
analyses in this paper to the LHC, but also more sophisticated searches using substructure 
techniques along the lines of [H] to look for boosted Z or h bosons. Also, the detectors have 
timing and pointing capabilities that should be exploited to their fullest to analyze the case 
of non-prompt NLSP decays [35] . 

( 7 )ln both ATLAS studies cited in this paper the gluino and squark masses for the model points studied 
were 0(1) TeV. 
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